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Abstract. A cosmological constant, A, is the most natural candidate to explain 
the origin of the dark energy (DE) component in the Universe. However, due to 
experimental evidence that the equation of state (EOS) of the DE could be evolving 
with timc/rcdshift (including the possibility that it might behave phantom- like near 
our time) has led theorists to emphasize that there might be a dynamical field (or some 
suitable combination of them) that could explain the behavior of the DE. While this 
is of course one possibility, here we show that there is no imperative need to invoke 
such dynamical fields and that a variable cosmological constant (including perhaps a 
variable Newton's constant too) may account in a natural way for all these features. 



1. Introduction 



The phenomenon of the accelerated expansion of the universe is presently one of the 
central issues of both observational and theoretical cosmology. A number of diverse 
cosmological observations [TJ |2] have by now established the accelerated nature of the 
present expansion of the universe and even provided additional information on the 
deceleration/acceleration transition and the redshift dependence of the expansion of 
the universe. From the theoretical side, the sole fact that the universe is presently 
accelerating, and may continue to do so, has triggered many studies. Some particularly 
interesting possibilities include braneworld models of the late-time cosmic acceleration 
The real theoretical challenge, however, lies in understanding the dynamics leading 
to the accelerated expansion of the universe. Despite the fact that many promising 
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models have been proposed, the fundamental nature of the accelerating mechanism 
remains presently unknown. The attempts towards shedding some light on the cause of 
the acceleration of the universe employ a broad range of concepts in many theoretical 
frameworks. The most widely used and the conceptually simplest option assumes the 
existence of dark energy (DE), the cosmic component with the negative pressure. Dark 
energy is a very useful concept since it encodes all our ignorance on the acceleration 
of the universe in a single cosmic component. Furthermore, DE can also be used as 
an effective description of other mechanisms of the acceleration of the universe [4 . A 
traditional candidate for the role of DE is the cosmological constant (CC), A. In fact, 
the cosmological FRW model with cold dark matter and A as the DE component (the so- 
called ACDM cosmology) fits the data reasonably well PjJ Theoretically, however, the 
CC as a DE candidate faces a formidable problem related to the very many (55 at least) 
orders of magnitude difference of its predicted value in quantum field theory (QFT) and 
the observed value [5] . This huge discrepancy leads indeed to the cosmological challenge 
of the millennium and calls for a more profound treatment of the CC problem jHj. The 
many inconsistencies related to the "A conundrum" have led to the development of the 
dynamical DE models. The tacit assumption of all these models is that A does not 
contribute to the DE or that it vanishes, which merely circumvents the CC problem. 
The dynamical DE models, however, incorporate the advantage that they approach the 
modeling of the mysterious dark energy component in a more general way, allowing 
its properties to vary with the expansion. Notwithstanding this welcome feature, the 
necessity of a (severe) fine-tuning of the parameters of these models in order that the 
measured value of the DE coincides with the predicted value does persist and it is in no 
way less worrisome than in the CC case. The list of dynamical DE models comprise, 
among others, quintessence [7] , phantom energy [H] , Chaplygin gas [§] etc E] . Very 
often these models of DE are realized in terms of dynamical scalar field(s). These fields 
were actually introduced on more or less phenomenological grounds to deal with the CC 
problem long ago [TQ1E]. In particular, in the popular quintessence approach the scalar 
fields and their parameters/scales are totally unrelated to the known particle physics 
fields, and as a consequence an obvious connection to fundamental physics is lacking. 
Not only so, in all these models the scalar field potentials just concocted to describe an 
acceptable form of DE are non-renormalizable by power counting and, what is worse, 
they are essentially field-theoretically unmotivated. 

In this paper we take the point of view that these DE models are not more 
fundamental than the original A model. However, we exploit the most useful property 
of the former, namely the dynamical character of the DE. Therefore, we generalize the 
CC concept allowing the variability of the A term (and possibly also of the gravitational 
coupling G) with the cosmic time, both assumptions being perfectly compatible with 
the cosmological principle and can be realized, as we shall see, in a fully covariant way. 
The support for such a generalization of the A term comes from QFT on curved space- 
time [121 EH| and/or quantum gravity approaches [Hj. Here, however, we do not derive 
the variability of A and G from these models, but discuss the general implications of 
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this variability. Our conclusions will therefore be valid for any specific model of the 
mentioned type. Let us emphasize that this approach embodies many virtues: i) the 
problem of the CC is dealt with instead of circumventing it; ii) the variability of A 
may shed some light on the observed value of A; iii) a variable A is indeed a type of 
dynamical DE, and will be handled here in the language of the effective DE picture. We 
will show in a precise way how a cosmology with variable parameters (A, G) can mimic 
a dynamical DE model and produce a non-trivial effective equation of state (EOS), 
Pd = w c sPd, beyond the naively expected one for the A term (w\ = —1). Since the 
EOS parametrization is widely used in all planned experiments aiming at a precise study 
of the DE [4| , such correspondence can be useful and it may even unveil some unexpected 
limitations of the EOS method for characterizing the ultimate nature of the DE. 



2. Dark energy picture versus variable cosmological parameters picture 

Before presenting the general procedure for obtaining the effective dark energy EOS 
corresponding to a generic variable A and G model, we briefly discuss the frameworks 
of two approaches which we call for short the DE picture and the variable CC picture. 
The DE picture assumes the existence of two separately conserved cosmological energy 
density components, the matter-radiation component and the DE component. In this 
picture the Einstein field equations 

Rnv - ^9^ R = 8itG Tfa, (1) 

have a total energy-momentum tensor of the form = T^ u + T®. The standard 
energy-momentum tensors of matter-radiation and the one for dark energy T® v are 
conserved separately. In the framework of FRW metric, the conservation of the energy- 
momentum tensor for matter-radiation, V M T* = 0, leads to the standard conservation 
law 

d ^ + aH DPs = 0. (2) 

Here a = 3(1 + u s ) where u s = and 1/3 for nonrelativistic matter and radiation, 
respectively. In the standard DE picture the assumption that V^T^ = results in the 
additional conservation law for the DE: 

^ + 3(l + w cS )H DPD = 0. (3) 

Generally the parameter u> eff is redshift dependent, u> e fj = w^(z). The Hubble 
parameter in this picture is defined by the Friedmann equation (assuming zero spatial 
curvature) 

Hh = ^(p s + PD ), (4) 

where a subscript D has been appended to H to distinguish the expression (J3J) from its 
counterpart in the variable CC picture - see below. The solution of the conservation 
laws (J2J) and (j2J) results in the following scaling laws for the components of the model: 

p s {z) = p s (0) (1 + zf (5) 
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p D (z) = p D (0)((z), (6) 




and 



l + z' J 

Using the expressions (0) and (|B|L the Friedmann equation acquires the form 

Hl{ z ) = Hl [n° M (i + *)« + n°, c (*)] • (7) 

On the other hand we have the variable CC picture, which describes the kind of models 
studied in this paper. Any of these models incorporates the matter-radiation component, 
a variable A component and, possibly although not necessarily, a variable gravitational 
coupling G. The variable A model represents a modification of Einstein equations 
of General Relativity which maintains its geometrical interpretation. The dynamical 
equations for gravity are given by 

R^ v ~ ^9fiuR = 8-itGT^ + g^A = SuGiT^ + g^ u p A ) , (8) 

where T^ u stands for the energy-momentum tensor of matter. This equation 
demonstrates that the full covariance is maintained even if A and G acquire space-time 
variability. Here p A is the energy density associated to A. In the framework of the FRW 
metric, p A and G depend on the cosmic time only, in accordance with the cosmological 
principle. Friedmann's equation is straightforwardly obtained from (JHJ) and reads 

Hl= 8 -^( P + PA ). (9) 

The general Bianchi identity of the Einstein tensor leads to the covariant conservation 
law 

VqG(T^ + WA )] = 0, (10) 
which for FRW metric acquires the form 

j t [G(p + p A )]+3GH A (p + p) = 0. (11) 

This "mixed" conservation law connects the variation of p A , G and p, where the scaling 
of the matter-radiation density p may be non-canonical in this picture (in contrast 
to p s in the DE picture). In this paper we consider a very broad class of models, 
just assuming the variability of the aforementioned quantities, without specifying the 
fundamental origin of such a variation. A number of variable CC models of various 
kinds [in] , and the renormalization group (RG) models of running p A and G [121 EH IH] 
provide the basis for the variability of these cosmological parameters. For example, the 
RG models [121 HB] do not determine the time dependence of p A and G directly, but 
indirectly specifying them in terms of other cosmic dynamic quantities (matter density 
p, Hubble parameter H, etc): 

Pk{z) = p A (p(z), H(z), ...) , G(z) = G(p(z), H{z), ...) . (12) 

These functions usually have a monotonic dependence when expressed as functions of 
cosmic time or redshift. The relations p2jl . with the general conservation law pijl . lead 
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to the complete solution of the variable CC cosmological model. Using this solution the 
general expression for the Hubble parameter becomes 

H\{z) = Hl [n M f M (z;r)(l + zr + n A f A (z-r)] . (13) 

Here fu and /a are known functions of redshift which may also depend on parameters 
r = ri,r 2 ,... originating from the fundamental dynamics. They generally have a 
nontrivial dependence on z and only in the case of ACDM cosmology they satisfy 
fu — /a = 1- Furthermore, they also fulfill the conditions /m(0, r) = 1 and /a(0, r) = 1 
in accordance with the cosmic sum rule Q° M + f2 A = 1- Notice that in general the two 
sets of cosmological parameters in the two pictures (|T3*j) and (J7J) will be different, e.g. 
AQm = &m — &m 7^ 0, because they correspond to two different fits of the same data. 

3. Matching of pictures and effective dark energy equation of state 

The two pictures presented in the preceding section may be considered as two separate, 
general DE models. In the remainder of the paper we, however, assume that they are 
equivalent descriptions of the same cosmological evolution. More precisely, we study 
the effective DE dynamics associated to the variable CC model through the procedure 
named the matching of pictures. The matching of pictures requires that the expansion 
history of the universe is the same in both pictures, i.e. that their Hubble functions 
are equal, Hp = at all times. In this way, for a known dynamics of the variable 
CC model, an effective DE can be constructed. A number of general results for the 
behavior of the effective DE density can be obtained with interesting implications to 
the observational data. The aforementioned matching of the two pictures gives the 
following equality connecting the dynamical cosmological quantities in both pictures, 
G(p + pa) = G (p s + po). Using H olt = —dzj (1 + z), the general Bianchi identity (fTTj) 
can be cast in the form 

(1 + z) d(p s + p D ) = a (p s + p D - £a) dz , (14) 

where we have introduced 

Uz) = ^Pk{z). (15) 

Next we use the scaling law (j2J) for p s and we arrive at the compact form for the redshift 
evolution of the effective DE density in the variable CC picture: 

d Pp( z ) _ a Pd( z ) ~ £a<» ( 16 ) 
dz 1 + z 

The integration of ()16|) readily yields a closed form expression for the effective DE 
density 



p D (z) = (l + zY 



p D (0) - a 

Jo 



o (l + ^') (a+1) 



17) 
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Expanding £a(#) around z = — 1 one can see that Pd(z) — > £a(z) at sufficiently late time 
(i.e. when z — ► —1). Finally, the effective EOS parameter for the variable CC model 
easily follows from (JHJ) and f)lfi|) 

«^,) = -1^(1 + ,)J^ = -1 + H(i-«£>), (18) 
3 p D dz 3 \ polz) / 

and it is seen to depend on the structure of the quantities (|T3|) and (|T7jl. From this 
expression we immediately see that in this kind of models w e s can cross the w e s = — 1 
line as soon as £a equals po- This observation reinforces the role of the effective DE 
density p^ in the study of the CC boundary crossing in the dark energy picture ||. 



4. Effective quintessence and phantom behavior of cosmological models 
with variable A and G 

For the following considerations it will be convenient to recast the solution of the ODE 
()lb|) in a different way. Let z* be the redshift value at which £a(^*) = Pd{ z *), i-e. 
w e s{z*) = — 1. Then it is easy to show that (fTtjj) is solved by the density function 

Pd{z)=U(z)-(1 + z) l {1 + z , )a dz> ■ (19) 

Quite remarkably, one can prove that a value z* always exists near our present time; 
namely in the recent past, in the immediate future or just at z* = 0. The proof of this 
assertion is obtained by straightforward calculation using: i) the matching condition 
Hy) = H\ and ii) the differential constraint that the general Bianchi identity 
imposes on the functions /a and fu in (|13j). In this way one obtains the following 
relation 

dp D (z) a(l + z) a ~ l 



dz 1 - Q° M 



-p D (0)(Q° M f M (z;r)-Q° M ) . (20) 



From this the proof is obvious. Indeed, we have already mentioned that at the present 
epoch function f M satisfies the condition f M (0,r) = 1. Moreover, the parameter 
difference AQm = — should not be large because the two pictures (DE and 
CC) describe the same physics. Hence from the continuity of fu it is clear that 
there must be a point z* close to where (1201) vanishes, that is to say a point that 
satisfies w c s(z*) = —1. The advantage of the formulation (fHJj) becomes evident when 
one calculates the slope of the po function: 

dpp(z) a _! r* dz' d^(z') 

— ; = —a (l + z) / — — — . (21) 

dz K ' Jz* (1 + z') a dz' K ' 

This compact expression reveals some counterintuitive and general aspects of the 
effective DE density evolution for variable CC models in which £,\(z) is a monotonous 
function of z. Intuitively one would expect that for £a(z) growing/ decreasing with 
z (decreasing/growing with expansion) pu should be quintessence-like/phantom-like. 

|| For other recent theoretical approaches to the w c g = — 1 boundary crossing see e.g. |18j . 
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The expression ()21|) shows that this is not the case. Namely, for £a(z) growing with 
z, pr> decreases with z for z > z*, i.e. in this redshift interval has the phantom-like 
characteristics. Only for z < z* po behaves as quintessence. Analogously, for £a(#) 
decreasing with z, pr> behaves like quintessence for z > z*, whereas only for z < z* it 
becomes phantom-like. Interestingly enough, if z* lies near our recent past (see Fig.QJi), 
the last case could just reflect the experimental indications of DE phantom-like behavior 
suggested by recent analyses of the data|20]. A concrete framework describing this last 
possibility is outlined in Section El These results illustrate that in variable CC models, 
the behavior of effective DE density is generally not determined by the CC only, but by 
the joint behavior of all quantities entering the general Bianchi identity pi)) . Especially 
interesting results are obtained when in the variable CC models the matter component 
p is separately conserved, i.e. when it also satisfies (J2J). In this case (fTTjl implies 
d^A/dt = —(p/Go) dG/dt, which from (|2Tj) results in 

= «(1 + ^ ^ P{z) - G(z*)\ . (22) 

Thus in this case the properties of pn depend only on the scaling of G with redshift, 
e.g. if G is asymptotically free and z* > z, then p D behaves effectively as quintessence 
(dpo/dz > 0), whereas if G is infrared-free then pr> behaves phantom-like (dpp/dz < 0). 

5. Effective dark energy picture of the RG model 

As a concrete illustration of the general procedure we have developed for obtaining 
the effective DE properties, in this section we consider the analysis [19) of the 
renormalization group cosmological model of ^B] characterized by G = const and 
the scaling law pa = C\ + C 2 H 2 for the CC. Here C 1 = pa,o — (^Hq)/(8ttG) and 
C*2 = (3z/) / (8irG) , where v is the single free parameter of the model -a typical value is 
\v\ = u = 1/127T [T9]. This model is fully analytically tractable and relatively simple 
expressions for u> eff can be obtained. In this particular case it is clear that (fT3j) reads 
£a(z) = Pa{z)- Therefore, for the flat universe case the effective parameter of EOS 
obtained from ()17J) and (JTSj) is 

„*(,) = -! + ( i-„)_ nS,(i + ^--'-^(i + ^ . (23) 

For \v\ 1 we may expand the previous result in first order in v. Assuming 
AO M = M° M - Q° M = we find 

w cS (z) ~ -l-3z/^|(l + ^) 3 ln(l + ^). (24) 

This result reflects the essential qualitative features of the general analysis presented in 
the previous sections. For v > 0, Eq. ()24j) clearly shows that we can get an (effective) 
phantom-like behavior (w c s < —1) and for v < we have (effective) quintessence 
behavior. We see that this variable CC model can give rise to two types of very different 
behaviors by just changing the sign of a single parameter. However, one can play with 
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Figure 1. (a) Numerical analysis of the effective EOS parameter w e s, Eq. (|23|l . as a 
function of the redshift for fixed v = — vq < 0, and for various values of Afijvf- The 
Universe is assumed to be spatially flat (f^ = 0) with the standard parameter choice 
fl° M = 0.3 , fl° A = 0.7; (b) Extended z range of the plot (a). We see that for Afl M < 
there exists a transition point z* near our recent past: namely, the one corresponding 
to the crossing of the CC barrier w e s — — 1 by the lowest curve in the figures. 



more parameters if desired. Indeed, as we have seen the cosmological parameters in the 
two pictures (DE versus CC picture) will generally be different (Afijf ^ 0). Fig[T] shows 
in a patent manner that in this case, even for v < 0, the variable CC model may exhibit 
phantom behavior due to the existence of a transition point z* in our recent past. 

6. Conclusions 

We have shown that a model with variable A (and possibly of G) generally leads 
to a non-trivial effective EOS, thus mimicking a canonical dynamical DE model (i.e. 
one with conserved DE density). That EOS can effectively appear as quintessence 
and even as phantom energy. Moreover, we have proven that there always exists a 
transition point z* near z = where w e s(z*) = — 1. If this point lies in our recent past 
(as illustrated in Fig.[IJi) there could have been a recent transition into an (effective) 
phantom regime w e s(z) < —1, as suggested by several analysis of the data. If, however, 
experiments would unambiguously indicate w e s(z) > —1 we could equally well interpret 
this in our framework, for z* could just be in our immediate future. We conclude that 
variable (p\, G) models may account for the observed evolution of the DE, without 
need of invoking any combination of fundamental quintessence and phantom fields. 
The eventual determination of an empirical EOS for the DE in the next generation of 
precision cosmology experiments should keep in mind this possibility. 
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